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Rate Dependence and Role of Disorder in Linearly Sheared 
Two-Dimensional Foams. 

Gijs Katgertjj Matthias E. Mobius, and Martin van Hecke 
Kamerlingh Onnes Lab, Universiteit Leiden, Postbus 9504, 2300 RA Leiden, The Netherlands 

(Dated: March 3, 2008) 

The shear flow of two dimensional foams is probed as a function of shear rate and disorder. 
Disordered, bidisperse foams exhibit strongly rate dependent velocity profiles. This behavior is 
captured quantitatively in a simple model based on the balance of the time-averaged drag forces in 
the system, which are found to exhibit power-law scaling with the foam velocity and strain rate. 
Disorder makes the scaling of the bulk drag forces different from that of the local inter-bubble drag 
forces, which we evidence by rheometrical measurements. In monodisperse, ordered foams, rate 
independent velocity profiles are found, which lends further credibility to this picture. 

PACS numbers: 47.57.Bc, 83.50.Rp, 83.80.Iz 



Similar to other disordered materials such as (colloidal) 
suspensions, granular media and emulsions, foams, which 
are dispersions of densely packedeas bubbles in liquid, 
exhibit a non-trivial rheology [l|, H H, 0, H, H, Q- When 
left unperturbed, foams jam into a metastable state cor- 
responding to a local minimum of the surface energy, 
where surface tension provides the restoring force under- 
lying their elastic response for small strains [g, |7|, [Sj . Un- 
der a continuous driving force, foam bubbles overcome 
these local minima and the foam starts to flow, and the 
viscous dissipation that arises in the thin fluid films that 
surround the gas bubbles becomes important. 

It is a daunting task to translate the bubble-bubble in- 
teractions to the rheology of foams [3|, |j, |9|, |lO| . Already 
for a single bubble sliding past a solid wall, Bretherton 
showed that the drag force scales nonlinearly with the 
bubble velocity [J, luJ, ll3] , and by analogy one would ex- 
pect the drag forces arising between sliding bubbles to 
be nonlinear also. In addition, at the multi bubble scale, 
foam flows can be disordered and intermittent [SJ, |8|, Isl . 
Foams share this combination of nonlinear interactions 
and complex flows with other disordered media. How- 
ever, bubble interactions are probably simpler than those 
of, e.g., frictional grains, and are similar to those of the 
soft spheres without static friction that have been stud- 
ied extensively in the context of jamming IJ, [ij, Il5l |. 
making foams eminently suited for fundamental studies 
of the flow of disordered media. 

In this Letter, we will address the role of disorder for 
foam flows, by experiments on the rheology of foams both 
at the coarse grained and at the bubble level. To probe 
and visualize foam flows, a number of experiments have 
been conducted recently in quasi two-dimensional geome- 
tries. Here the foam flow is driven by moving sidewalls, 
and the soap bubbles either form a bubble raft where they 
freely float on the fluid phase [l6|j_are sandwiched by two 
glass plates in a Hele-Shaw cell 17|, or are trapped be- 
tween the fluid phase and a top-plate [l8|, ,12,] . The pres- 
ence of such a top-plate leads to shear banding of the 
flow hjl. This can be understood from the additional 




FIG. 1: (a) Schematic top view of the experimental setup, 
showing how two counter rotating wheels partially immersed 
in the fiuid and spaced by a gap W shear the foam. Data 
is taken in the highlighted area and a typical fiow profile is 
indicated, (b) Side view showing the layer of bubbles trapped 
below the top plate and the grooved shearing wheels, vq is 
the a;-component of the wheels angular velocity, and is equal 
to (jjro over the contact line (dashed) , since vq — cjri cos (j> — 
■ ujro. The applied strain rate -ya equals 2vo/W . 



drag forces exerted on the bubbles flowing under the top 
plate, which will be balanced by gradients in the bulk 
stresses of the material. 

A model based on the balance of drag forces which 
captures the observed shear banding qualitatively was 
recently introduced by Janiaud et al. [IQ] . For simplicity, 
it was assumed that the drag forces exerted by the top 
plate scale linearly with bubble velocity, and that the 
bulk stress varies linearly with strain rate. These linear 
laws lead to rate independent flows [lO|. 

Here we experimentally probe the flow of such 2D 
foams which are trapped between the fluid phase and 
a top-plate. We find that the flow depends crucially on 
the applied strain rate ja- disordered, bidisperse foams 
exhibit rate dependent flow profiles, which become in- 
creasingly shear-banded for large 7^. 

These findings are captured in a model in which the 
time-averaged drag forces between bubble and top plate. 




■ V(,= 0.026 mm/s. 

■ Vo= 0.26 m m/s - 
Vp=2.6mm/s 
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FIG. 2: (Color online), (a) Experimental image of bidisperse foam (left) and corresponding bubble tracks (right). Note the 
swirling motion, (b) Rescaled velocity profiles v/vo for W = 7 cm and no as indicated, compared to profiles obtained from our 
model Eq. Owith a = 2/3, /3 = 0.36 and k = 3.75 (red curves). Inset: v/vq for vo — 8.3 mm/s and W equal to 5, 7 and 9 cm 
— for convenience, we chose the origin at the left boundary here, (c) Illustration of the model defined by Eq. ([5]). 



Ffc^, and between neighboring bubbles, Fbi, are balanced. 
While the continuum limit of our model is similar in spirit 
to the model of Janiaud, the crucial new ingredient are 
non-linear scaling laws for the wall drag and the bulk 
stress — these non-linear scalings are essential for cap- 
turing the observed rate dependence. 

We establish the precise scaling forms of the averaged 
drag forces in disordered foams by varying the applied 
shear rate over three orders in magnitude and fitting the 
data to our model, and confirm these scalings by indepen- 
dent rheological measurements. We furthermore perform 
rheometrical measurements on ordered lanes of bubbles, 
which reflect the viscous drag force between individual 
bubbles. Surprisingly, the averaged drag forces in the dis- 
ordered foam scales differently from the local drag forces 
between individual bubbles, and in our range of param- 
eters, the averaged forces are much larger than expected 
from naively scaling up the local drag forces. In contrast, 
for monodisperse, ordered foams, the local drag forces, 
averaged drag forces and top-plate drag all scale simi- 
larly, causing rate- independent flows [l9|. We attribute 
the modiflcation of the drag forces to the disordered and 
non-affine motion of the bubbles in the bidisperse foam 

Setup — A bidisperse (50:50 number ratio) bubble 
monolayer is produced by flowing nitrogen through two 
syringe needles immersed at flxed depth in a soapy solu- 
tion consisting of 5 % volume fraction Dawn dishwashing 
liquid and 15 % glycerol in demineralized water (viscos- 
ity 77 = 1.8 ± 0.1 mPa-s and surface tension cr = 28 ± 1 
mN/m). The resulting bubbles of 1.8 ± 0.1 and 2.7 ± 0.1 
mm diameter are gently mixed to produce a disordered 
bidisperse monolayer and are covered with a glass plate 
(see Fig. [T|). The weighted average bubble diameter (d) 
is 2.25 mm. 

Two parallel PMMA wheels of 195 mm radius and 
9 mm thickness are partially immersed in the liquid 
through 10 mm wide slits in the top plate such that 



they are in contact with the foam over a length of 230 
mm, while having an adjustable gap distance W ranging 
from 50 to 100 mm (Fig. 1). The wheels have a rough- 
ness of order 3 mm at the contact line due to etched 
grooves, like the spokes on a bicycle wheel, to ensure no 
slip boundaries for the bubbles, and are counter-rotated 
by two micro-stepper motors. The bubbles bridge to the 
top plate, which is completely wetted by the soap solu- 
tion, and we fix the liquid fraction of the foam by keeping 
the distance between glass plate and liquid surface fixed 
at 2.25 ± 0.01 mm. Coalescence, segregation and coars- 
ening as well as the drag force between bubbles and fiuid 
phase are negligible. 

The average velocity v{y) in the x-direction is ob- 
tained from both particle tracking and particle image 
velocimetry-like techniques. Since the time-resolved flow 
is strongly disordered and intermittent (Fig. 2a), we av- 
erage over time and over x, where we restrict the x-range 
to a central region of length 60 mm (Fig. la) where re- 
circulation is neghgible {{vy) — 0). 

Rate dependent flows — We measured the flow profiles 
V for gap width W equal to 5, 7 and 9 cm, and driving 
velocities va = 0.026, 0.083, 0.26, 0.83, 2.6 and 8.3 mm/s. 
In Fig. [2] we show a few examples of these. The main 
observation is that the velocity profiles strongly vary with 
the driving velocity vq, and become increasingly shear 
banded for large vq. 

Drag force balance model — The flow profiles and the 
scaling forms of the drag forces are connected by a sim- 
ple model in which the average drag forces are balanced. 
As illustrated in Fig. [2j;, we divide the foam in lanes of 
width {d) and balance the time-averaged top plate drag 
per bubble F^,^ with the time-averaged viscous drag per 
bubble due to the lane to the left (Ff,;,) 
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FIG. 3: (a) _F(,„ is deduced by trapping a monolayer of bub- 
bles between a rough bottom and a smooth top plate (in- 
set). From the power-law scaling of torque T we deduce the 
drag force per bubble as a function of Ca [Jl, and find that 
Fi„ = /6„Ca°-'5^±°°2, with f^^ ^ i.5±0.1xl0-3 N. (b) Ft^ is 
deduced from the time-averaged torque exerted on our bidis- 
perse foam as a function of ACa{= rjAv/a) (filled circles). 
The foam is sheared in a Couette cell of inner radius 1.25 cm, 
outer radius 2.5 cm (hence a gap of 5 bubble diameters) with- 
out a top plate (inset). We obtain Fbb ~ fY+.fbb{ACa)" , with 
the yield threshold /y « 1.2(5) x 10"^ N,fbb ~ 5.6(9) x 10"'' N 
and /3 = 0.40(2)) (solid line). Open circles are the same data 
with the yield torque obtained from the fit subtracted, which 
are well fit by a pure power-law with exponent 0.4 (dashed 
line), (c) Drag force between two pinned and ordered bubble 
lanes in a Taylor-Couette geometry. The black line indicates 
power-law scaling with exponent 2/3. 



strongly, we assume that we can express the average drag 
forces in terms of the average velocities ^'^ We non- 
dimensionalize velocities according to the definition of 
the capillary number {Ca :— rjv/a), and propose: 



F 



bw 



Fbb — 
^ bb — 



fbUvv'/<^r , (2) 

fY + fbb[{ri/<j){v'-v^-'f , (3) 

fY + .fbb[{v/<^W^'-ylf ■ (4) 



The expression for Fbw is essentially the result for a sin- 
gle bubble sliding past a solid wall, for which Bretherton 
showed that the drag force Fi,^ scales non-linearly with 
the capillary number 



m 12,^ 



fbw is a constant 
with dimensions of force of order arc, where Vc is the 
radius of the bubble- wall contact [20| . The power-law in- 
dex a depends on the surfactant. Dawn has a low surface 
shear modulus [2l|, for which a = 2/3 [3] (see Fig. 3a). 

For Fbb we conjecture a Hcrschcl-Bulkley type expres- 
sion, which combines a finite threshold fy with a power- 
law dissipative term. The crucial exponent (3 will be 
determined from the flow profiles and rhcology below. 

Inserting these expressions into Eq. ([T]) and defining 
k = fbw/ fbb we arrive at: 
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where it should be noted that the yield threshold fy 
drops out of the equations of motion — we keep it here 



to remain consistent with our rheological measurements 
(see Fig. 3). 

Model vs. experimental flow profiles — To compare our 
model (Eq. [5]) to the eighteen experimental flow profiles 
obtained for three widths and six driving velocities, we 
need to determine the two dimensionless parameters /? 
and k. To avoid being affected by edge effects near the 
shearing wheels, we focus on the part of the data where 
\v\ < 3/4 -uo, and solve Eq. ([5|) by numerically integrating 
from where f = to the y value for which v = 3/4 • vq. 
For fixed P and k we can thus compare the experimental 
data and model prediction. 

To determine /3 and fc, we require that all profiles 
are fit well for the same values of these fitting param- 
eters. When /3 is not chosen optimally, we find that k 
systematically varies with vq, but for (3 — 0.36 ± 0.05, 
this systematic variation is minimized. We find that for 
a = 0.67, /3 = 0.36 and k = 3.75 ± 0.5, all 18 data sets 
can be fit excellently by our model (Fig. 2) [22l |. 

Constitutive Relation — By taking the continuum limit 
we can rewrite our model Eq. ([S]) as the balance between 
a ID body force and gradient of the ID stress r: 



2/3 



dr 



^ / n . \ 0.36 

r = Ty + fbb ' 



(6) 

(7) 



where Ty is an undetermined yield stress. This yields the 
constitutive equation for a Herschel-Bulkley fluid, and 
the value /3 — 0.36 is remarkably close to recent results 
for 3D bulk rheology of emulsions and foams [2, |j| . 

Rheological determination of a, (3 and k — The force 
laws that underly our model can be probed directly 
by rheological measurements, and we have measured 
the bubble-wall and inter-bubble forces with an Anton 
Paar MCR-501 rheometer (See Fig. 3a-b). We find that 
a = 0.67 ± 0.02, thus confirming that Fbw is given by 
Bretherton's law for mobile surfactants. The measured 
value of the exponent (3, 0.40 ± 0.02 is within error bars 
to what we found by simply fitting the model to the fiow 
profiles. We extract from the rheological measurements 
an estimate for the ratio k = fbw/ fbb ~ 2.5 ± 0.5. This is 
close to the value k = 3.75 ± 0.5 estimated from the flow 
proflles ,23). 

Rheological determination of (3 at the bubble scale — 
To see if the effective inter-bubble drag force corresponds 
trivially to the drag forces at the bubble scale, we have 
probed this force by narrowing the gap width of our Cou- 
ette geometry so that two pinned and perfectly ordered 
lanes of bubbles slide past each other, and find that here 
Tbb ~ (Aw)2/3 (Fig. [2;)(2i|. We believe that this result 
refiects the actual viscous drag force between individual 
bubbles sliding past one other, which implies that the av- 
erage drag forces in a bidisperse, disordered foam derive 
in a highly non-trivial way from the drag forces at the 




FIG. 4: (a) Experimental image of monodisperse foam (left) 
and bubble tracks (right). Note the order in the foam and the 
absence of swirly motion in the tracks, (b) Velocity profiles 
for a monodisperse foam (d — 2.7mm) at 7 cm gap, for 0.083 
mm/s (black), 0.26 mm/s (dark grey) and 0.83 mm/s (light 
grey). Red curves are fits to the model with both a and /? 
equal to 2/3, and k = 0.3. 



bubble level. 

Discussion — The drag forces exerted on the bubbles 
by the top plate, which at first sight might be seen as ob- 
scuring the bulk rheology of the foam, enable us to back 
out the effective inter-bubble drag forces and constitutive 
relation of foams from the average velocity profiles. By 
comparing these results with rheometrical measurements, 
we note a remarkable difference between the scaling of the 
drag forces at the bubble level and the bulk level: we find 



(Av) 



0.36 



at 



Fbb ~ (AvY'-^ at the bubble level and Ftt 
the bulk level. 

One might understand this anomalous scaling as fol- 
lows: The degree of disorder does not affect the drag 
forces at the bubble scale, but it does modify the bub- 
ble motion. For disordered foams, the bubbles exhibit 
non-affine and irregular motion — hence they 'rub" their 
neighboring bubbles much more than when they would 
flow orderly, and consequently, the averaged viscous dis- 
sipation is enhanced over what could naively be expected 
from the local drag forces [9|]. This picture is corrob- 
orated by recent simulations on the bubble model [3|, 
where one recovers this "renormalization" of the drag 
force exponent 25|, |26[ and rate-dependent flow profiles 

To further illustrate our picture, we linearly shear a 
monodisperse foam (d = 2.7 mm) and recover rate inde- 
pendent profiles (see Fig. 2]). By tracking we confirm the 
absence of signiflcant disordered bubble motion in this 
case. Our model only yields rate independent profiles if 
a — (3 so that (3 = 2/3, which implies that without disor- 
der, the connection between local and bulk drag forces is 
trivial. This also solves the conundrum why Wang et al. 
[19| found rate- independent fiow profiles, as their foams 
are essentially monodisperse. 

In conclusion, polydisperse, disordered foams exhibit 
rate dependent fiows due to anomalous scaling of the 
averaged drag forces Fbb- We suggest that anomalous 
scaling of bulk properties caused by non-affine motion at 
the particle scale may be a general feature of disordered 



systems close to jamming [3|, |9|, [isl, [ij, [iS | . 
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